Magnetic disorder and dynamical properties of a Bose-Einstein condensate in atomic 

waveguides 
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We systematically investigate the properties of the quenched disorder potential in an atomic 
waveguide, and study its effects to the dynamics of condensate in the strong disorder region. We 
show that even very small wire shape fluctuations can cause strong disorder potential along the wire 
direction, leading to the fragmentation phenomena as the condensate is close to the wire surface. 
The generic disorder potential is Gaussian correlated random potential with vanishing correlations 
in both short and long wavelength limits and with a strong correlation weight at a finite length scale, 
set by the atom-wire distance. When the condensate is fragmentized, we investigate the coherent 
and incoherent dynamics of the condensate, and demonstrate that it can undergo a crossover from a 
coherent condensate to an insulating Bose-glass phase in strong disorder (or low density) regime. Our 
numerical results obtained within the meanfield approximation are semi-quantitatively consistent 
with the experimental results. 



I. INTRODUCTION 

Low dimensional physics has long been an important 
and extensively studied subjects in condensed matter 
physics since the early theoretical interest in 60th. In 
addition to the ordinary solid state system, e.g. semi- 
conductor quantum wells, quantum wires, carbon nan- 
otubes, and organic conductors etc., systems of ultracold 
atoms in highly anisotropic magneto-optical traps have 
been also a promising new system for the low dimensional 
physics of both fermion and boson particles^*^. Among 
varies proposals of the anisotropic magneto-optical traps, 
microfabricated magnetic trap, also called atomic waveg- 
uide or microtrap^, has an additional important advan- 
tage for the coherent transport of atoms along the quasi- 
one-dimensional (Q1D) waveguide potential (for a typ- 
ical system setup, see Fig. QJa) and the description in 
the next section). Such combination of ultracold atoms 
in quantum mechanical limit with the great versatility 
of the fabrication technique has opened a new direc- 
tion to study many Q1D quantum physics, like coherent 
transport, beam splitter—, interference of matter wave£, 
and Tonks gas limitS. Many theoretical works have been 
proposed in the literature to study these aspects in very 
recent years. 

From experimental point of view, however, in order to 
successfully trap ultracold atoms in the atomic waveg- 
uide, one has to apply very strong confinement potential 
in the transverse dimension, which is usually achieved by 
moving atoms very close to the conducting wire on the 
substrate. (Electric current in the wire cannot be too 
large in order to avoid heating 9 .) When the atom- wire 
distance is smaller than some critical values, which is usu- 
ally about 100—200 fim, however, nontrivial atom density 
modulation (fragmentation) of the condensate occurs, as 
have been observed by different groups&i2iiiii2ii^. It is 
also found that atom could cannot be transported with- 
out excitations to higher energy band when such frag- 
mentation occurs inside the clouc&ii. For a static frag- 
mentized atom cloud, it shows a rather generic (but not 



universal) length scale, A, in z direction, which becomes 
smaller when d is reduced. Kraft et. al^ further observe 
that the positions of atom density maximum/minimum 
can be exchanged if the direction of the offset magnetic 
field is reversed, showing a nontrivial quenched disorder 
effects in the microtrap system. Such interesting disorder 
effects can be very crucial when considering the coherent 
transport in the magnetic waveguide, and was first inves- 
tigated theoretically by us recentlji 2 ^. We note that such 
quenched disorder induced fragmentation phenomena is 
different from the thermal fluctuation studied in Ref. 26 . 
The thermal fluctuation effects can be neglected in our 
present scope of interest, because it becomes prominent 
only when the atoms are even closer to the wire surface 
(say d < 20 fJxa), while the fragmentation phenomena 
occur at well-larger distance (say d < 200 /im). 

The disorder effects in an interacting bosonic system 
have been an interesting subject in the context of liq- 
uid helium system for decades. In weak disorder limit, 
it is believed that the disorder effects is irrelevant to the 
ground state properties due to the strong repulsive in- 
teraction between bosons in two and three dimensional 
systemsi£ii£iiL±&. Similar conclusion also applies to one- 
dimensional (ID) bosonic system 19 . However, when the 
disorder strength increases, it is argued by Fisher et al. 
that the system can undergo a quantum phase transition 
to a "Bose-glass" phase 20 , which breaks the usual U(l) 
symmetry of bosons to have a condensate but is insu- 
lating due to disorders. However, the existence of the 
"Bose-glass" phase in the liquid helium system is still 
unclear in the present experimental data2i. 

However, the invention of optical lattice for trapping 
ultracold atoms provides another route to study the dis- 
order effects to the bosonic systems. Experimentalists 
uses laser of different frequency to produce an artificial 
random potential (speckle pattern) and study the ground 
state or transport properties of cold atoms22i2i. Recently, 
the onset of Bose glass regime is observed near the Mott- 
insulator regime^i. 

In this paper, we extend our earlier work 2 ^ and pro- 
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vide a more detailed study on the origin, properties, and 
effects of disorder potential in atomic waveguide system. 
We show a first principle and quantitative theory to de- 
scribe the disorder potential in the atomic waveguide, 
and demonstrate that even a very small wire shape fluc- 
tuation can cause a strong disorder field, which has a 
correlation function that vanishes at small and infinite 
wavevectors and is peaked at a wavevector, k c with a 
length scale determined by the atom- wire distance (d). 
We note that similar analytical works and numerical com- 
parison with experimental data have been also studied 
recently2L2& after our work. In strong disorder regime, 
the condensate becomes fragmentized as observed in the 
experimenljSiifliiiiiSiii. We then concentrate on the dy- 
namical response of the condensate in the shaking exper- 
iment and propose that a quantum phase transition from 
a superfluid state to an insulating Bose glass state may 
be observed in the parameter regime of current experi- 
ments. Our results provide a good starting point of the 
disorder effects to the coherent transport and to the low 
dimensional Bose systems in the future development. 

This paper is organised as follows: In Section [n] we 
present the theory of the magnetic disorder potential gen- 
erated from the wire shape distortion in the microchip 
system. In Section IIIII we study the condensate frag- 
mentation in strong disorder region. Then we discuss 
the finite temperature effects and other details in Sec. 
IIVI and summarize our results in Section 



II. DISORDERED MAGNETIC FIELD 

In Fig. ^a) we show a typical experimental setup of 
an atomic waveguide: ultra-cold atoms are loaded into 
the microtrap, which is composed by a radiant magnetic 
field gradient (generated by the electric current in a one- 
dimensional microfabricated copper wire) and a uniform 
bias field, B±, in the transverse direction (y). Atoms 
are confined at the potential minimum with a distance 
d = to the wire center (in c.g.s. unit), where B± 
exactly cancels the field of the wire current. A uniform 
offset field, B\\, is applied parallel in the wire direction 
(z) to reduce the trap loss by polarizing the atom spins. 
Additional magnetic field gradients can be also applied 
in z direction to confine the atoms in the longitudinal 
direction. The elongated confinement potential can be 
approximated by a harmonic potential, with confinement 
frequencies, uj± and uu in x — y and z directions respec- 
tively. 



A. General properties of disordered magnetic field 

Before studying the disorder magnetic field generated 
by the wire shape fluctuations, we think it is helpful to 
study some more general properties of the disorder mag- 
netic field. Using the fact that the typical Zeeman energy 
of atoms is very large (~ MHz), we can safely assume that 



all atoms are condensed in the lowest spin state (i.e. fully 
polarized in the direction of the total magnetic field) so 
that the confinement potential, U (r) , is proportional to 
the total magnetic field: 

C/(r)=/i |B tot (r)|, (1) 

where /i is the magnetic dipole moment of atoms and 
B tot (r) =B\\z + B ± y + B (r) + 5B(r); SB(r) is the dis- 
order magnetic field, which will be calculated in details 
later. Since the condensate in the microchip system is 
highly elongated in z direction with a very small trans- 
verse radius, Rtf ~ 1 — 5 /im, we can neglect the finite 
size effects for simplicity and just consider the disorder 
potential at the confinement center (x,y,z) — (d, 0, z), 
where Bj^ = B±y cancels the unperturbed azimuthal 
field, Bo- In the limit of small disorder magnetic field, 
one can make an expansion of Eq. JQ| and obtain: 

XTTf\ , RMj ^a[\mz)\ 2 -25B z {zf] 
6U(z) ~ s\\n a SB z (z) + ,(2) 

where sy ± 1 is the sign of the offset field, B», in z di- 
rection. Note that in Eq. J3J) the first (dominant) term 
is linearly proportional to SB Z , so that SU(z) changes 
sign if the direction of either the electric current or the 
offset magnetic field is changed. This simple observation 
explains the experimental results observed in Ref. |12| . 
where Kraft et al. find that the positions of local poten- 
tial maximum/minimum can be exchanged by changing 
the direction of By or the current. The second term of 
Eq. (J2J, however, does not change sign for different By or 
current direction, and therefore explain why the density 
profile of the condensates obtained by opposite directions 
of current (or opposite B») are not symmetric. However, 
since \Bu | is in general much larger than the disorder 
field (|5B/2?||| <C 0.1% in general), we will simply ne- 
glect the second term of Eq. J3J) and consider SB Z only 
when applying to the realistic disorder field calculation, 
although we will still derive the disorder magnetic field 
in all the components of <5B later, which may be useful 
in considering the finite size effects of BEC in the future 
study. 

In this paper, we consider a microchip system where 
the magnetic field is generated from a rectangular cooper 
wire with width Wo (in the substrate plane) and hight 
H (vertical to the substrate plane, see Fig. Ha)). The 
typical scale for the semiconductor wire is Ho ~ 1 — 5 
fxm, and Wo ~ 3 — 50 /xm^ii. In general, the wire shape 
fluctuation can occur in both of these two directions due 
to the fluctuations during the crystal growth. However, 
in this paper we assume the width fluctuation or center 
changes only in the horizontal direction (y), but not in 
the vertical direction (x) for simplicity (see Fig. [TJ. This 
is a good approximation because the current density in x 
direction should generate much smaller magnetic field at 
the position right above the wire according to the Biot- 
Savart law in classical electrodynamics. Such approxi- 
mation is also confirmed by the numerical work in Ref. 
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FIG. 1: (a) Schematic figure of a BEC loaded in a micro- 
trap. Bo(r), -B||£, and B±y are the azimuthal magnetic field 
generated by a steady current, the offset field, and the bias 
field respectively, (b) and (c) present the C— and W — types 
of wire shape fluctuations respectively. 



[28j after our first analytical studyS^. Therefore we can 
define the wire shape function, y = —Wo/ 2 + Jl{z) for 
the left boundary and y = Wq/ '2 + fn{z) for the right 
boundary of the wire in the substrate {y — z) plane, and 
the current fluctuation becomes in y — z directions only, 
i. e. 



SJ(y,z) ~ 5J y (y,z)y + SJ z (y,z)z, 



(3) 



which are closely related to the shape of the wire and will 
be studied in details below. 

When studying the current fluctuation in the wire, it is 
in general tentative to approximate the disorder nature of 
current fluctuation by the following expectation values^ 



(SJ(r)) dh 





5 a pD(r ■ 



(4) 
(5) 



where a, (3 = x, y, z, (■ ■ -)dis is disorder ensemble average, 
and D(r — r') is the correlation function of the fluctuat- 
ing current. However, we point out that Eq. JSJ is not 
correct for a normal conduction wire, because it assumes 
no correlation between different components of the fluc- 
tuating current. For any static current distribution, as 
we will see below, the microscopic continuity equation 
still applies to the disorder current so that any current 
fluctuation in one direction should somehow affects the 
current in other directions. Besides, the total current 
following in any single wire (i.e. for any specific system 
as done in the realistic experiment) must be a constant 
too, not only the ensemble-averaged value is unchanged 
as implied by Eq. Therefore the realistic current 

fluctuation should obey some more strict conditions than 
those shown in Eqs. I0J and JSJ. This fact will change 



the disorder nature of the magnetic field significantly as 
we will show in the rest of this section. 



B. Current conservation and the boundary 
conditions 

The first step is to introduce the current conservation 
and the electro-dynamical equations in a self-consistent 
method. We start from the static version of continuity 
equation, V • SJ(r) = 0, which gives 



dSJ y (y,z) d5J z (y,z) 



dy 



dz 



= o. 



(6) 



Another equation can be derived from the Maxwell equa- 
tion, V x E = -dB/dt, with dB/dt = for static cur- 
rent and assuming a constant conductivity, ct, inside the 
wire (and zero outside the wire) (i.e. J = crE), so that 
V x J = or equivalently: 



dSJ z {y,z) dSJ y (y,z) 



dy 



dz 



0. 



(7) 



Eqs. 10 and J7J) are the two equations we have to solve 
by incorporating the proper boundary conditions. 

The first boundary condition is from the total current 
conservation in z direction, which gives (See Fig. E c )) 

pHo/2 r-W /2+f B (z) 

Io = dx' dy'3(y',z)-z 

J H /2 J-W /2+f L (z) 

= H / dy'(J + 6J z (y' : zj) 

J-W a /2+f L (z) 

~ Io + H J (f R (z)-f L (z)) 

r-W /2 

+H / dy'SJ z (y',z), (8) 

J-Wo/2 

where Joz is the unperturbed (or average) current density 
in the wire, Iq — JqWqHq is the total current, and in 
the last equation we have expand the right hand side by 
assuming small wire shape fluctuation, i.e. 



\h/R(z)\^W 

In other words, we have 

-W /2 



(9) 



Jo (/«(*)- = - 



Wo/2 



dy'SJ z (y',z). (10) 



The second boundary condition can be obtained by re- 
quiring the current flow is parallel to the wire bound- 
ary, whose tilted angle from the z axis is given by 
dfL/w( z )/dz for the left and right boundary. Therefore 
we have (see Fig. IHb)) 

df L / R (z) = SJ v (TW /2 + f L/R (z),z) 
dz Jo + SJ z (TW /2 + f L/R (z),z) 

^ SJ v ( T Wo/2,z) ^ 
Jo 
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where we have used the approximation of small wire 
shape fluctuation (Eq. ©) and assumed the shape fluc- 
tuation is smooth, i.e. 



\fL/ R i*)\ « 



1. 



(12) 



(Of course we always have |5Jy, a | <C Jo-) It is easy to 
show that Eq. I|ll[l can directly imply Eq. I|10(l if we take 
the z derivative of the later and use the continuity equa- 
tion, Eq. ©. This shows that Eqs. ©,0, and Jill) are 
the three independent equations we have to use to solve 
this problem self-consistently. The current conservation 
and boundary conditions has been included fully in the 
limit of small wire shape fluctuation (Eq. and l|12|)). 



C. Solution of current density fluctuations 

To solve above equations, Eqs. ©-(01 we define an 
auxiliary function Gj(y',z') which gives the fluctuating 
current density as follows 



5J v (y,z) = J ( 



dGj{y,z) 
dz 



(13) 



dGj(y,z) 
dy 



(14) 



Therefore the current continuity equation, Eq. JBJ, has 
been satisfy automatically, and then Gj(y,z) must sat- 
isfy a Laplacian equation, V 2 Gj(y, z) = inside the wire 
according to Eq. 0. We can use the standard method 
of separation of variables: 



fe>o 



x (Cfc cos(fcz) + D k sin(fcz)) , (15) 
so that the boundary conditions (Eq. (fTTjO become: 



Gj(±W Q /2,z) = f R/L (z). 



(16) 



Defining f R/z (z) = £ £ fe>0 f*J, L cos(kz) + f^ L sin{kz) 

we can solve Gj(y, z) via Eq. I|16(l after a straightforward 
algebra: 



cos(fcz) 



isinh(Wo) 



sinh(W /2) cosh(fc?/)/ 1 c fc + cosh(fcWo/2) sinh(fcy)/- 



sin(fcz) 
sinh(fcVFo) 



smh(kW /2) cosh(fcy)/ 2 c fc + cosh(fcW /2) smh(ky)ff k 



(17) 



where f^ k = |(/^ 2)jfc ± /f (2);fc ) is the cosine(sine) 
Fourier components of the wire center/ width fluctuations 
(see Fig. 0b)-(c)). Therefore we can solve the charge 
density fluctuation directly by taking the derivative as 
shown in Eq. (|14f) in the limit of small wire fluctuations. 



obtain the disorder vector potential, <5A(r), to be 



D. Disordered magnetic field 

Applying the approximation of small wire shape fluc- 
tuation and above results for current fluctuations, we can 



J 



5A(r) 



A(r)-A (r) 

I fHo/2 



dx > r dz > r l2+tR{z) iy ,i^mmmm - Mr 

' •/-Ho/2 J-oo J-W /2+f L (z) y/(x - x') 2 + (y - y') 2 + (z - z') 2 

H '' /2 dx > r dz > r 2 iy> , *UMv + 6U M z 

' J-H /2 J-oo J-W /2 yj{x ~ X 1 ) 2 + (y - y') 2 + (z - Z 1 ) 2 
Ho/2 



1 



1 



poo 

dx' I dz' 



c J -H /2 



f R (z)J z 



!l{z)J z 



yj(x - x') 2 + {y- Wo/2) 2 + (z - z') 2 yflji: - x') 2 + (y + W /2) 2 + {z - z' f 



5 



Jo 



*> r dz > r /2 jy> g >w>*w*-w-<»-™ 

c J-H/2 J-oo 



W0/2 {{x - x') 2 + (y- y') 2 + {z- z') 2 ) 3/2 ' 



(18) 



where we have used Eq. (|16fl to cancel the extra terms via integration by parts. The magnetic field can then be 
obtained directly: 



<SB(r) 



00 J , f Ho/2 J / r Wo/2 J ,n t * >A 2 ( x - x 'f - (v - y'f - ( z - z 'f . 

dz' dx' dy'Gj(y',z')< — — - V^i 

C ./. .™ / TJ ~ In I t,t/_ In t t nc — <y> l\2 _1_ <„, — ti'\2 1 I ~ _ ~l\2\ t >/ Jj 



Jo 



-00 J-Ho/2 J-Wo/2 

-3(x-x')(y-y>) 



((s - x') 2 + (y- y') 2 + {z- z') 2 )°' 
— 3(a; — x')(z — z') 



((x - x 1 ) 2 + {y- y') 2 + (z- z') 2 f /2 {{x - x') 2 + (y - y') 2 + (z - z') 2 f' 2 ' 



r 



(19) 



Note that Eqs. @, lfT7)l and JT^J are the main results 
of this section, using the following approximation: (1) 
the wire boundary fluctuation in the vertical (x) direc- 
tion can be neglected, and (2) the wire shape fluctua- 
tions in y directions are small, i.e. \Jl/r( z )\ <S Wo an d 
\f' L / R (z)\ *C 1. We believe these are reasonable approxi- 
mations in the regime of present experimental interest. 



further approximate it by using the fact that (i) the wire 
hight, Hq, is usually much smaller than d in the x di- 
rection, and (ii) in practice we just need to study SB Z 
at (x,y,z) = (d, 0, z) according to the earlier discussion 
about Eq. (J2J). We then evaluate 8B z (d, 0, z) by integrat- 
ing the wire width and obtain 



E. Finite wire size effect and disorder correlation 
function 

Although one can numerically ca lculate the disorder 
field from Eq. I|19(l as done in Ref. [27l28j . but we may 



6B z (d,0,z) 



tt 7 r+°° 



, f w ^ 2 , -3Gj(y',z')d(z-z') 
c J-00 Z J-Wo/2 V (d 2 +y' 2 + (z-z') 2 Y/ 2 
6H J 



cd 



dz' 



(z-z') 



E 



(-l)"5-7--(2n + 3) 



(1 + (z - z') 2 f/ 2 ^ nl ■ 2«(1 + [z- z') 2 ) n J Q 



W /2d 



Gj{y'd,z'd)y' 2n dy', (20) 



where we have changed to a dimensionless dummy variable inside the integration. The last integral can be evaluated 
directly to be (the odd part of y in Gj(y, z) has been integrated out to be zero) 

"""w** - \ e w/& + *q*m T(2 " + - kw " /2) - 7(2 " + h ^ m 



sinh(Wo) 



2{kd) 



2n+l 



(21) 



Jq x n 1 e x dx is the incomplete 



where 7(71, x) ' ••" 1 - 

Gamma function. Therefore we can calculate the 
disorder correlation function, (5U(z)SU(z'))dis, whose 
Fourier component in momentum space, = 
/ dz e-' lk ^- z 'YSU{z)5U{z')) dis , is 



where 



D(x,y) 



2 S inh(a;/2) y> (-1)" K n+1 (y) 
xsinh(a:) n\ ■ (2y) n K x {y) 

x [ 7 (2n + 1, x/2) - 7 (2n + 1, -x/2)] (23) 



27oMa\ (kd) 4 



d 4 



\D{kW Q ,k d )Ki{kd)\ 2 F k ,{22) 



and F k = fdze- ik ^-^{f c (z)fc(z')). K t (x) is the 
modified Bessel function of the second kind. 
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FIG. 2: (a) Calculated disorder correlation function, A*,. 
Solid and dashed lines are for ffe=constant and Fj, oc e n 
with rj = d/2 respectively. For comparison, we also show a 
Gaussian correlation function (dotted line) of length scale d. 
(b) Typical disorder potentials in real space, calculated from 
the disorder correlation function shown in (a). Note that the 
scales of the vertical axes in both figures are not the same for 
different curves, up to the overall strengths of disorders. 



F. Disorder strength 



To measure the strength of disorder, we define the 
strength u s to be u 2 s = (u(z)u(z)) = J |§A fc , which is 
obtained by integrating the whole spectrum of the corre- 
lation function, A&. We then obtain 




1 



dk{kd) 4 'K 1 {kd) 2 D{kWQ,kdfF k 



1 f°° 

J dpp^KiipfD 

S(W ,d;{F k }) 
nd 3 



( pWp 
{ d ■ 



F, 



p/d 
(24) 



where p = kd is dimensionless variable, and a new func- 
tion, S(Wq, d; {Fk}), is defined to incorporate the wire 
shape disorder, which is proportional to F in large d 
limit. Therefore, if considering Wo <C d limit, we may 
simplify above result further to be 



= C- 



cd 



J) 1/2 = c». B± . (§)'%, 



where C is of order of one and depends on Wq and dis- 
order length scale only weakly for large d. According to 



u : d = 50 fim 

o d=100 

: d=150 

6 1 d = 200 

4 ■'"'/ 

2 L 



2 3 
k X (100 itm) 




FIG. 3: (a) Calculated disorder correlation function, Aj., for 
different atom-wire distance d. (b) Calculated density profiles 
of a 23 Na condensate by using the same parameters as the ex- 
periments (Fig. 4 of Refi-). Dotted lines are the density with- 
out disorder potential. In both figures, Fk is approximated 
by the form described in the text. 



Eq. i|25fl. for a fixed current (Jo), the disorder strength u s 
scales as d~ 2 5 for large d, roughly consistent with the ex- 
perimental result shown in Rcf. 1 -, where they estimated 
the disorder field ~ d~ 2,2 by equating the chemical po- 
tential to the disorder field at the onset of the fragmenta- 
tion. Besides, Eq. (|24|l shows that even if the wire center 
fluctuation is very small (Fq <C d 3 ), it can still generates 
a large disorder potential compared to the BEC chemi- 
cal potential, /i, because the energy of the bias magnetic 
field, /J, a B±, is in general of order of MHz, much larger 
than /i, which is just of order of kHz. Therefore it seems 
almost impossible to reduce the disorder effect in the mi- 
crochip experiment by improving the wire sample, and 
one has to consider the disorder effect to the coherent 
transport more seriously especially when d is not large. 



G. Model wire fluctuation and numerical results 

In Fig. OJa) we show the numerical results of A^ 
with Wq = for simplicity (the results are similar for 
Wq < O.ld). Solid(dashed) line is for Fk =constant(Jfc oc 
e _fc "' 7 with rj/d — 0.5). Results of a disorder potential 
with a Gaussian correlation function (i.e. A^ oc e _fc d , 
which cannot be a result of a current conserving calcula- 
tion as we discussed earlier) is also shown in the same 
figure (dotted line) for comparison. Unlike the usual 
Gaussian correlation most adapted in the literature, A^, 
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obtained in Eq. I|22[). is zero at fc = and peaked at 
a finite wavevector fc ~ 1.33/d (solid line) even if the 
wire center fluctuation has zero correlation length (i.e. 
Ffc=constant). In other words such length scale can be 
completely generic and is mainly determined by the mi- 
crotrap geometry rather than by the length scale of the 
underlying surface disorder of the wire. This interesting 
results are due to the current conservation inside the wire, 
which requires the current fluctuation to be cancelled in 
low wavelength limit (kd -C 1). The finite width of the 
peak shows that such disorder can be considered as a 
Gaussian correlated random potential. In Fig. EJb), we 
show typical disorder potentials in real space based on 
the disorder correlation functions discussed above. One 
can see that the usual Gaussian-type correlation (dotted 
line) cannot give any periodicity in the disorder poten- 
tial. On the other hand, if we assume a Gaussian corre- 

.2 2 

lation for the wire center fluctuation, i.e. Fk oc e n 
(dashed lines), it will help to reduce the contribution 
from wavevectors higher than 1 /d, and makes the quasi- 
periodicity more transparent as shown in Fig. |2Ib). The 
peak position of is then shifted to lower value fci, 
which gives the length scale, A = 2n/ki, of the Gaussian 
correlated random potential in real space. 

Another interesting situation is that if the wire shape 
has not only one random disorder but also an intrin- 
sic periodic disorder with a period r/o (i.e. Fk ~ 

Ci5k,k + C2e~ h n , where fco = ^/Vo and Ci/ 2 are rel- 
ative strengths of these two disorders), Eq. H22|) then 
has a "double-peak" structure in A& at k = fci and hi- 
The coexistence of these two length scales (no, which is d- 
independent and A, which is d-dependent) in a BEC frag- 
mentation can explain the double periodicity observed in 
Ref. When temperature is raised above T c , the 

chemical potential may be smaller than one (stronger) 
disorder but still larger than the other (weaker) one, so 
that only one length is revealed in the fragmentation of 
a thermal cloudii. 

Although the exact form of the wire shape fluctu- 
ation (fc(zj) has to be determined by a microscope 
as have been done in Refs. and is sample- 

dependent, we think it can be generally described by 
assuming F k = ^( e -^- k ^ 2 ^ + e-( fc + fc i) 2 ^), where 
(fc( z )fc( z )) = 2rnVn measures the wn " e shape fluctu- 
ations. For k\ — 0, we have a normal (Gaussian) dis- 
tribution function of fc( z ) with a length scale 771. For 
kii]i > 1, fc(z) has a Gaussian correlated structure with 
a length scale, A = 2n/ki. For 771 = 0, we recover the zero 
length scale fluctuation (white noise) for fc{z) and then 
k\ = ko ~ 1.33/rf as shown in Fig. |2Ja). In this paper, 
we choose our parameters to model the fragmentation 
observed in MIT group£ and will set (fc(z)fc(z)}={0.1 
/Ltm) 2 , rji = 100 /jm, and 27r/fci ~ 200 /jm for the numer- 
ical calculation in the rest of this paper. 

In Fig. Ota) we show the numerical results of A*, for 
different d by using above model fluctuation of the wire 
shape. As expected, the underlying disorder (now has a 



length scale in the Gaussian correlated form) will change 
the position of the peak of A&, while the overall length 
scale still depends on d. In Fig. B^b) we calculate a static 
BEC density profile of sodium atoms in different disor- 
der potential strengths by solving the ID static Gross- 
Pitaevskii equation (GPE): 

Q = _dg* + [ v+m+ |$|2_ M ]$ (26) 
2m 

within Thomas-Fermi (TF) approximation 2 ^, where m is 
the atom mass, and Vq(z) = ^muj^z 2 is the axial har- 
monic confinement potential. 

din = -pr, (aj_ = J— — ) (27) 

is the effective ID interaction strength incorporating the 
radial confinement potential 8 '^: a is the 3D s-wave scat- 
tering length ujj_ is the single particle confinement poten- 
tial in the transverse mode and a± is the typical radius 
in the transverse direction. We can see that the con- 
densate becomes fragmentized when d becomes small, as 
observed in the experiments. 

III. FRAGMENTATION IN STRONG 
DISORDER LIMIT 

In this section, we focus on the strong disorder regime, 
where the condensate is fragmentized into several pieces 
with weak coupling strength between the neighboring 
fragments We will use the modified tight-binding model 
to study the dynamical properties of these condensate 
fragments in the superfluid phase regime by using the 
dynamical mcanfield (Gross-Pittaivskii) theory. 

A. Tight-binding approximation 

For the fragments observed in the atomic waveguide 
experiment aVW 2 , it seems quiet reasonable to assume 
that the single particle tunneling between two neighbor- 
ing fragments is so weak that the system can be described 
by a tight-binding approximation. However, unlike the 
usual tight-binding model used in the optical lattice, the 
atom-atom interaction effects is very crucial in the frag- 
ment case in determining the onsite particle wavefunc- 
tion. We start from the full second quantization repre- 
sentation of ID boson problem: 

i-^ = -^ + (V a + SU- f i)4, + g 1D ^^, (28) 

where i[>(z, t) and V>(z, t)^ are bosonic operators. Now 
we introduce the onsite single particle wavefunction to 
describe the discrete nature and expand the single boson 
operator as follows 

$(z,t) - («)&,•(<), (29) 

3 



8 



where cf>j(z) is the onsite groundstate wavefunction of 
the local potential well j, which is centered at z = Zj. 
a,j is the bosonic operator for that specific eigenstate. In 
Eq. (|29[) we have assumed that all the boson atoms in 
each well are in the local groundstate, and no higher en- 
ergy modes can be excited (single band approximation) 
in each well. We will discuss the validity of such approxi- 
mation later. Substituting above equation into Eq. I28|) . 
and integrating the z coordinate after multiplying (j>j(z) 
in both sides of the Eq. I)28|l , we obtain (after neglecting 
the next nearest hopping energy and interaction between 
nearest neighboring wells) 



da-j 
'1H 



- ^2 K j,j+a a j+a + UjO^Clj 



a = ±l 



dz 



a,j, 
(30) 



where /\ ,. , = -2 J dzd z (f>j(z)d z <f>j>(z) + J dz<f>j(z) 
tfiji(z)(Vo(z) + SU(z) ~ fi)] is tunneling amplitude be- 
tween jth to j'th well. Uj = <7id / rf^|<^j| 4 is the onsite 
charging energy. 

Note that if we use the noninteracting ground state, 
(j)jo, for each well to expand the single particle operator 
in Eq. (|29|) . we obtain the well-known single-band Bose 
Hubbard model (BHM>2i: 



da-j 



■ 2^ Kj t j +a aj +a + (V ,j 
a=±l 



SUj 



+U j a]a J a j 



fj,)aj 
(31) 



where Vqj = J dzVo\(f>jo\ 2 and SUj = J dzSU\cf>jo\ 2 is the 
onsite smooth potential and disorder potential. (The lo- 
cal kinetic energy, i.e. the first term in the second line of 
Eq. i|3T|) . can be absorbed into the chemical potential.) 
The approximation of such noninteracting wavefunction- 
based single band BHM is self-consistent only if the onsite 
excitation energy, w; oc , is much larger than any other en- 
ergy scale in Eq. l|3*TJ) . Kj,j+i, SUj, Vbj Uj, and NjUj, 
where Nj = (aja^) is the average onsite number of atoms. 
The last criterion (u>i oc 3> NjUj) is required to confirm 
that the atom-atom interaction does not change the on- 
site wavefunction from the noninteracting one, 4>j0j which 
determines the values of Kjj+i, Uj, and SUj etc. as 
shown above. If uii oc is of the same order as NjUj, the 
interaction between atoms in each well can deform their 
wavefunction from the noninteracting ground state, V'jCb 
by including higher energy onsite eigenstates. As a result, 
the single band approximation used in Eq. (|29|l fails. 

However, in our present condensate fragment system, 
the number of particles in a fragment is so large (Nj ~ 
10 4-5 ) that the criteria, u>i oc NjUj cannot be satisfied 
in most of the experimental situations (for the fragment 
size, we can calculate that typically Uj ~ 0.05 Hz and 
uioc ~ 10 — 50 Hz). As a result, it is still reasonable to 
assume that atoms in each fragment can be still described 



by a single (condensate) wavefunction, cfij, which, differ- 
ent from the noninteracting wavefunction </>jo, is deter- 
mined mainly by the competition between the atom-atom 
interaction (instead of the local kinetic energy) and local 
confinement provided by the disorder potential. There- 
fore we can assume the U(l) symmetry is broken in each 



well, so that the onsite wavefunction, 



can be deter- 



mined by the following local static GP equation: 

= - 2^V 2 0, + (V +5U- [i)tfij + 3lD iVjV,| 2 ^,(32) 

where <pj is normalized to one and N® is of the equilib- 
rium number of particles in well j. Using above equation 
to eliminate the last term of Eq. Ij30|l , we obtain 



da j 
'~dt 



Kj,j+ a dj+a 

a=±l 



(33) 



Note Eq. is different from Eq. lETT)) . because the 
onsite wavefunction now have to be determined by in- 
cluding interaction, local kinetic energy, and local trap- 
ping potential via Eq. Ij32(l . As a result, the criterion 
of single band approximation used in the regular BHM 
(Eq. lO) can be softened to be lui oc S> Kjj+x,Uj for 
Eq. (|33JI . At the same time we still assume the temper- 
ature is low enough that the whole system is in quasi- 
condensate region, i.e. local density fluctuation in each 
well and the phase difference between neighboring wells 
can be neglected in determining the equilibrium onsite 
wavefunction. The randomness of the disorder potential 
is now absorbed into the the randomness of N® and also 
<j)j, which makes Uj and Kjji are also random. We note 
that a simplified Hamiltonian of strong disorder in both 
diagonal and off-diagonal term (i.e. without interaction 
Uj of Eq. (|3*3l ~) has been also studied recently in Ref. 

m 

The dynamics of Eq. can be studied within 

its meanfield version by taking (&j(t)) = (a]-(i)) = 

Nj(t)e lSj ^\ where Nj and Sj are the number of par- 
ticle and phase function of the jth condensate fragment 
respectively: 



N; 



E KjJ+WNify+aMSj+a - Sj) (34) 



a=±l 



9 K id+a\l 
1 a=±l V 

-(Nj-N^Uj 
-(Nj-N^Uj, 



Na 



cos(S j+a - Sj) 



(35) 



where we have neglected the single particle tunneling 
term for simplicity, which can be shown to be much 
smaller than the onsite charging energy (see below). It 
is interesting to note that Eqs. 134|) - (|35|l are very simi- 
lar to the well-known equations for Josephson junction^. 
However, the interaction between bosonic atoms leads to 



9 



an effective onsite charging energy (the last term of Eq. 

which closes the equations by changing the num- 
ber of particles and phase functions simultaneously. Note 
that throughout this paper we will neglect the dynamics 
of the quantum fluctuation completely for simplicity, be- 
cause it is very small due to the large average number of 
particles in each fragment (Nf > 10 4 ). 



B. Probability distribution of parameters in the 
tight-binding model 

In Eqs. the parameters, Kjj + i, Uj and N°, 

are all random numbers, because the disorder poten- 
tial in the microchip may generate fragments of different 
sizes and at different (random) positions, which affects 
the tunneling amplitude and onsite charging energy as 
well. In our calculation, we use a Gaussian trial wave- 
function to solve the onsite wavefunction, 4>j{z), of Eq. 
(|32|l at each potential j, keeping the chemical potential 
fixed by the total number of particles of the whole con- 
densate. The single particle tunneling amplitude Kj t j+i, 
onsite interaction energy, Uj, and number of atoms in 
jth fragment N®, are calculated by <f>j(z) as in the stan- 
dard method^i. In principle these three quantities are 
not independent of each other, but it is still useful and in- 
structive to study their probability distribution function 
individually for a given disordered magnetic field studied 
in the previous section. In Fig. we show the distribu- 
tion function of them for relatively weak (s = 0.5 dotted 
line) and strong (s = 1.2 solid line) disorder potential, 
where s = u s /u s and u s is the disorder strength defined 
in Eq. I|24fl with the same disorder realization as used in 
Fig. at d — 100 /zm. This is calculated by considering 
2 x 10 6 number of sodium atoms in ten potential wells 
for each time and then averaging the results over different 
(more than twenty thousands) disorder realizations. 

Several features can be observed in Fig. [3j (i) when 
disorder is enhanced, the distribution functions of all pa- 
rameters become broadened as expected and also change 
their mean values, (ii) For the number of particles per 
site (Fig. Ufa)), stronger potential wells can trap more 
atoms, leaving only very small number (< 10 2 per site) 
of atoms in the weaker potential wells. Therefore in 
the fragment situation we are considering (Eqs. i|32|) - 
(|30[l b some weaker potential wells can have only very few 
atoms, which behave like a weak link between their two 
neighboring wells. This can be observed in Fig. Efb) that 
(iii) when the disorder strength increases, more "weak" 
junctions between fragments appear so that the whole 
system may becomes more closer to an insulating phase. 
Simultaneously, we find that (iv) the onsite charing en- 
ergy is also increased by disorder as shown in Fig. Ufc). 
This is due to the fact that stronger onsite disorder po- 
tential can reduce the size of the fragment, leading to 
stronger onsite interaction energy, (v) Finally, in Fig. 
Eld), we show the distribution function of the "effective 
onsite disorder potential", N®Uj, as derived in Eq. (|33fl . 
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FIG. 4: (a)-(d) are respectively the probability distribution 
of local number of atoms N? , tunneling energy, Kjj+i, onsite 
charging energy, Uj, and onsite effective disorder potential 
N°Uj (see Eq. fSJJ. Dotted (solid) lines are for disorder 
strength s — u s /u a s — 0.5(1.2) respectively. The units of 
energy variables are normalized to Hz (h — 1). 



The double peak structure in weak disorder limit (dot- 
ted line) is due to the different contributions of N° and 
Uj respectively. When disorder increases, they merge to 
one peak and move to stronger side. Therefore we can 
summarize that the disorder potential of the microchip 
fragmentation can increase the inhomogeneity of the den- 
sity profile, reduce the tunneling amplitude, Kjj+i, in- 
crease the onsite charging energy, Uj, and also increase 
onsite potential variation. All of these features lead to 
a strong indication of an insulating phase in strong dis- 
order region. In the rest of this section, We will discuss 
the dynamics of such system under the shaking experi- 
ment and study how it can be related to the superfluid 
to insulator transition. 



C. Two- fragment dynamics and self-trapping 

We start with the two- fragment case (j = 1,2), where 
some unique dynamical properties due to the nonlinear 
GPE can be described more precisely. For the two frag- 
ment case, the only eigen frequency is luj = V 2KU 
(K ee K ia y/N?N$ and U = \{U l + U 2 )) for the small os- 
cillation of AN = N 2 -N 1 ~ (iV 2 ° - ATf ) and AS = S 2 - Si . 
It is easy to show that Eqs. I|34|) - (|35|) are equivalent 
to the problem of single planar pendulum, where AS(t) 
and AN(t) represent the angular position of the weight 
from the vertical line and the angula r velocity respec- 
tively (we assume \fN\N 2 ~ y N^N® for small density 
variation). In such simple two- fragment situation, an in- 
teresting phenomenon, "self-trapping" effect can be ob- 
served in the two-fragment case when AS is larger than 
7T. As shown in the phase portrait of Fig. E{a), the sys- 
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FIG. 5: Dynamics of two identical fragment system: (a) Flow 
diagram in the phase space, (b) Time evolution of the density 
variation with different initial AN. (c) Oscillation frequency 
and (d) amplitude (about its average value) as a function of 
initial density modulation from equilibrium value. 

tern trajectories of the initial conditions denoted by the 
first two right arrows do not pass the origin (AN = 0) 
and keep flow away by growing AS exponentially. In 
Fig. Efb), we show a numerical results for AN(t) from 
different initial AiV(0) (with AS(Q) = 0), using a set of 
typical parameters: N® = N§ = 10 5 , Ki,2 = 1 Hz and 
U = 0.05 Hz. We can see that when initial displacement, 
AN(0)/N%, is less than 1% (solid line), the number of 
atoms oscillates between these two fragments with fre- 
quency luj = 100 Hz. When AN(0) increases, the oscil- 
lation amplitude first increases accordingly (dotted line) 
and then suddenly decreases for AN(Q) > AN C ~ 0.04iVi 
(dashed and dash-dotted lines). Besides, the density im- 
balance AN(t) then never passes zero but just oscillates 
about a new average value. Such counterintuitive result 
is due to the nonlinear nature of Eq. (|34|l . and can be 
realized in a single pendulum problem where the weight 
has a very large initial velocity to overcome the gravita- 
tion potential even at the highest position of the circle 
(AS = 7r) and then keep rotating forward with a nonzero 
velocity (see the trajectory following the first two right 
arrow in Fig. EJa)). It is also very similar to the Ac 
Josephson effect, where a constant electric potential drop 
between the two connecting superconductors can cause a 
oscillating current between them. In our present situa- 
tion, the initial potential drop is provided by the large 
imbalance of number of particles so that the fast oscilla- 
tion of single particle tunneling cannot reduce such initial 



imbalance except for other damping mechanism. We can 
calculate the critical amplitude of the initial density vari- 
ation easily and obtain 

\AN C \=2J-^. (36) 

We note that such critical density variation becomes 
smaller when the tunneling energy is weaker and/or the 
average onsite charging energy is larger, i.e. close to the 
insulating phase. 

In Fig. GJc)-(d), we plot the oscillation frequency and 
amplitude as a function of different initial density vari- 
ation, £(0) = AN(0)/N avc . It is easy to see that the 
oscillation amplitude drops very fast when initial dis- 
placement is larger than the self-trapping point. This 
result also exists even when considering the full quantum 
mechanics in such simple two well systemM, because the 
typical number of particles per site in the fragments is 
so large (> 10 4 ) that the quantum fluctuation can be 
neglected. 

D. Shaking experiment in condensate fragments 

In optical lattice system, the dynamics of a condensate 
cloud can be studied by the "shaking experiment" , a sud- 
denly shifting of the global confinement potential with 
a finite displacement, and then observing the successive 
center of mass motioniS^. When the initial displace 
is small, the condensate oscillates harmonically, showing 
a coherent Josephson junction tunneling between neigh- 
boring wells. When the displacement is larger than some 
critical value, however, the center of mass motion be- 
comes strongly damped, indicating a dynamical instabil- 
ity, which is a classical phase transition due to the break- 
down of meanfield solution 35 . It is believed that when the 
strength of optical lattice is tuned to be strong enough, 
the shaking experiment with small displacement can be 
used to investigate the proposed superfluid to (Mott) in- 
sulator transitions^, where (unlike in the superfluid 
phase) the small displacement should not result in a co- 
herent center of mass motion 3 ^. In our earlier work 25 , 
we proposed that similar experiments can be done to 
investigate the superfluid to insulator transition in the 
microchip fragments, where the insulating phase is best 
understood as a Bose glass phase due to the underlying 
disorder nature^i. Here we will study such multi-well 
dynamics in more details via the meanfield equation of 
motion derived in Eqs. l|33 jl -l|35 )l . 

However, although there are many similarity between 
the condensate fragments we consider here and the con- 
densates in optical lattice^^ 2 ^!, their difference in the 
sizes of local potential wells does bring some significant 
difference of their dynamics. We first note that due to 
the large size of the disorder potential well in the mi- 
crochip, the overall condensate density profile is far away 
from the result of Thomas-Fermi inverse paraba (com- 
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FIG. 6: Schematic potential and density profile of condensate 
fragments in a microchip system, (a) Self-trapping effects 
for large displacement, D. The atoms in fragment 3 cannot 
tunneling as a condensate into the potential well 4, which is 
generated by the large displacement of global confinement po- 
tential. The noisy curve in the upper portion denotes the local 
disorder potential, and the lower portion shows the fragment 
density profile before shaking. Solid/dashed lines denote the 
global confinement potential before/after shaking, (b) Den- 
sity variation for small displacement (no disorder potential is 
shown here) .The dashed line for the density profile denotes 
the new equilibrium profile in the new potential after shak- 
ing. Solid lines are the density profile before shaking, (c) 
The potential and density variation with respect to the new 
potential and equilibrium density profile. Such ANj becomes 
the source of fragment dynamics as time evolves. Up-down 
arrows indicate the possible density oscillation in each well. 



pared the solid and dotted lines in Fig. Efb)). Actu- 
ally, one can expect that the density profile has a kind of 
"discontinuity" around the edge of the whole condensate. 
Applying the self-trapping dynamics discussed in previ- 
ous section, we then expect that the atoms in the edge 
fragments will not effective tunneling into the potential 
well next to it, which is created by shifting the confine- 
ment potential (see Fig. Eta))- Such strong edge effects is 
right due to the nonlinearity-induced self-trapping phe- 
nomena as discussed earlier. Therefore we do not expect 
that the condensate fragments will be driven to motion 
as a whole by such shaking experiment. The only pos- 
sible exception is by thermalizing the atoms in the edge 
fragments to higher energy modes, which is certainly not 
a coherent motion at all and is not our current interest 
in this paper. 



Despite of such huge difference, the shaking experiment 
with small displacement can be still applied to study the 
coherent motion in the condensate fragments by realizing 
the density profiles of the whole condensate. As schemat- 
ically shown in Fig. Etb); although the small displace- 
ment D may not change the whole condensate position 
with respect to the disorder potential, it is indeed capa- 
ble to change the onsite potential strength to a slightly 
new value, and hence the atoms in each fragments be- 
come to flow between neighboring wells to response the 
change of local chemical potential. Precisely speaking, 
such small displacement shift should change the shape 
of local potential well also, and hence change their po- 
sition as well as well width/depth. However, we believe 
it is a reasonable approximation to neglect the change of 
well positions and widths, and concentrate on the effects 
mainly from the deviation of local chemical potential, 
which gives a new equilibrium density profile, A" etu , so 

that the old density profile, N° ld , becomes an initial non- 
equilibrium source for the coming dynamics. In Fig. EI C )- 
we schematically depict the linear change of (parabolic) 
confinement potential and hence the unbalanced density 
source for the future dynamics, ANj = N° ld - Aj ietu . 
If the whole condensate is in superfluid regime, we ex- 
pect that ANj(t) = Nj(t) - Nf ew will oscillate about 
zero with a definitely phase relation with respect to its 
neighboring site (vertical dashed arrows in Fig. Et c ))i 
when \ANj(t)\ is small enough. In general, the density 
variation {Q(t) = \Nj(t) - N$\/N$, where N° = N™ ew 
and Nj(t = 0) = N° ld ) is largest at the edge fragments, 
because the potential deviation is the largest there while 
its average number of atoms (A^ ) is in general smaller 
than the fragments in the center of condensate. 

In the following calculation of fragment dynamics, we 
shift the constant change of chemical potential so that 
AV Q {z) = f uj 2 {z + D) 2 - fuj 2 z 2 = muj 2 Dz = Vz is zero 
at the center of the whole condensate with an effective 
potential displacement T>. The initial density variation 
(i.e. the density before shaking, N° ld , respect to the new 
equilibrium density, N™ ew , after shaking) can be approx- 
imated by ANj(t = 0) ~ AVo(Zj)/Uj. We will therefore 
consider how the density variation (ANj(t)) evolves as 
a function of time with different displacement, T>. Be- 
sides, we also make a further approximation by assuming 
the global confinement potential can be neglected when 
calculating the onsite wavefunction and other dynami- 
cal variables, because they are mainly determined by the 
local disorder potential. 

E. Numerical results 

As an example of the condensate fragment dynamics, 
we show in Fig. [7|the typical results for a four fragment 
system. We choose K\ % — ^3,4 = 1 Hz ^2,3 = 0.5 
Hz, J7i, 2 ,3 = 0.05 Hz, and A? = A° = 1.3 x'lO 5 and 
= A3 = 1.2 x 10 5 as some typical values shown in 
Section llllBI Several features can be observed from Fig. 
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FIG. 7: (a) Density and (b) phase oscillation in each frag- 
ments of a typical four fragment condensate. In both figures, 
solid, dotted, and dashed lines are results for O=0.15, 0.2, 
and 0.3 Hz/(j,m respectively. In (a) the four sets of curves 
are for well j = 1, 2, 3, and 4 respectively from bottom to 
top. In (b) the three sets of curves are phase gradients in the 
neighboring sites: Sn — Si, S3 — S2, and S4 — S3 respectively 
from bottom to top. 



(i) When the displacement is small (solid lines), both 
density profiles and phase gradients oscillate about their 
new equilibrium values with small amplitude, (ii) When 
the displacement is increased, one of the phase gradient 
(here it is S3 — S%) becomes unbound, leading to a self- 
trapping phenomena where the density profile oscillates 
around a nonzero mean value (dashed lines), (iii) Be- 
tween the generic oscillation in small displacement and 
the self-trapping in the large displacement, we find that 
the condensate dynamics becomes chaotic in the interme- 
diate range of displacement (dotted lines). Such chaotic 
dynamics can be investigated via the frequency spectrum 
(see Ref. [2^) or the time correlation function. This is a 
classical instability of the meanfield GPE, and its result 
to the center of mass motion in optical lattice has been 
investigated in Refs. [35J. 

In Fig. [S] we plot a numerical calculated dynami- 
cal phase diagram for ten-fragment system with N tot = 
2 x 10 6 sodium atoms in the disorder potential calculated 
earlier (at d = 100 /im) in terms of disorder strength (s) 
and critical displacement, D. Here the dashed line sep- 
arates the multi-mode oscillation (e.g. solid lines in Fig. 

01 from the chaotic motion (e.g. dotted lines in Fig. |7J, 
while the solid line separates the multi-mode/chaotic mo- 
tion from the self-trapping motion (e.g. dashed lines in 
Fig. 0. We can see that the critical displacement of 
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FIG. 8: Critical shaking displacement of a ten- fragment sys- 
tem as a function of disorder strength (s = u a /u° a ). Total 
number of atoms is 2 x 10 6 . 



self-trapping decreases dramatically for s = u s /u s < 1.1 
and then becomes smoothly for s > 1.1. It never reaches 
zero in our meanfield approximation. We believe that 
after including the full quantum fluctuation, the mean- 
field solution becomes incorrect in the strong disorder 
region, where the quantum fluctuation can decrease the 
critical displacement so much that even infinitesmall dis- 
placement to the condensate will be self-trapped without 
coherent motion (i.e. the solid line of Fig. [SJterminates at 
finite disorder strength, say s = 1.1 This can be un- 
derstand as a signature of quantum phase transition from 
superfluid phase in small disorder to an insulating phase 
in strong disorder, which is best understood as a Bose 
glass phased due to the nature of randomness. However, 
since the low energy excitation properties in such strong 
disorder regime are still poorly understood, we could not 
exclude the possibility of different quantum phases, such 
as Mott glass as suggested by Giamarchi et. al~ due to 
the competition between a (white noise) random poten- 
tial and a commensurate periodic potential. 



F. Estimate quantum fluctuation effects 

The transition from superfluid phase to insulator phase 
is driven by the quantum fluctuation in the competi- 
tion between interaction and random potential 20 . In- 
cluding the quantum fluctuations in the dynamics of 
condensates has been studied in some limited cases«2isi 
by either solving the full coupled nonlinear Gross- 
Pitaevskii-Bogoliubov-de-Gennie equations ", using dy- 
namical meanfield variational method 31 , or by using trun- 
cated Wigner approximation^. Studies including ran- 
dom potential is much more difficult and still in progress. 
Here we will give some estimate about such transition us- 
ing the statistical properties, which, at least in principle, 
should be able to be probed by the dynamical shaking 
experiment when displacement is tuned to infinitesmall. 

Transition into the insulating state may be character- 
ized by the ratio of the "Josephson energy" between the 
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FIG. 9: The dotted line shows the probability P(Q < 1). 
The system is in the superfluid state when this probability is 
close to zero (see text), while it is in an insulating state when 
it is close to one. The solid line shows the probability that 
one can observe fluctuations in the number of atoms between 
the neighboring wells with fluctuations in the number of par- 
ticles larger than £ m i n = 0.1 and oscillation frequency larger 
than u m i n = 2ir x 1 Hz. When this probabilities approaches 
zero, the system will appear self-trapped in the shaking ex- 
periments. Here iV avc is the average number of atoms in a 
single well (mini-condensate). Disorder strength for (a) are 
the same as used in Fig. |3J) for d = 100 /im (denoted to be 
u l l), while it is five times stronger in (b). The results are av- 
eraged over more than twenty thousand pairs of two-fragment 
system with average number of atoms, iV ave , per well. 



neighboring wells to the charging energy of one of the 
wells, i.e. Qj = (jj'+u-lxyV • Without disorder poten- 
tial, the meanfield calculation 4 ^^ estimates the super- 
fluid (SF) to Mott insulator (MI) transition at Q = 1. 
For the Gaussian correlated random potential we discuss 
here, we expect that quantum fluctuation induced SF to 
BG phase transition 44 should also appear at roughly the 
same place. Due to the randomness nature of the conden- 
sate fragments, the effects of quantum fluctuation is im- 
portant if the probability to have Qj < 1 (=V(Q < 1)), 
is of the order of one, where the probability if obtained 
by averaging over disorder ensemble. In Fig. [5] we show 
the calculated V(Q < 1) as a function of the atom den- 
sity for two different strengths of the disorder potential 
(dotted lines). In order to capture the quantum fluc- 
tuation nature of each random junction, we calculate Q 
for each pair of fragments across the potential barrier 
with average number of atoms JVa Ve per well. According 
to Fig. El we can estimate that the quantum fluctua- 
tions are crucially important to a pair of fragments when 
iV avc < 10 4 - 4 - 2.5 x 10 4 for s = u s /u a s = 1. This critical 
average number of particles become higher (iV avc < 10 D ) 
for stronger disorder potential (s = u s /vP s = 5). In 
the same figure, we also show the calculated probability 



(solid lines) to have Josephson frequencies u>j to be larger 
than some minimum frequency, w m in = 2wx 1 Hz, and the 
density contrast, £ = to be larger than (,„;„ = 10%. 
The values of u> m m and £ m ; n are given by the experimen- 
tal resolution of frequency and density deviations^ for the 
density oscillation. We can see that this probability (de- 
noted by V(u> > cjmin, C > Cmin)) decreases when average 
number of atoms per well decreases and/or the disorder 
strength increases (toward the insulating phase). Both of 
these two results (V(Q < 1) and V{uj > u) m i n , £ > Cmin)) 
suggest that one can observe the superfluid to insula- 
tor (Bose glass) transition 44 in the parameter regime of 
present experiments. 

It is interesting to compare above estimate of quan- 
tum fluctuations (Fig. |§J) with the dynamical phase di- 
agram associate with the shaking experiments shown in 
Fig. 03 In Fig. |H| the dynamical properties are calcu- 
lated for an ensemble of ten fragments, instead of two. 
Since the superfluid dynamics of the whole condensate 
can be strongly suppressed if one of these junctions (say, 
junction I between fragment I and 1 + 1) becomes insulat- 
ing due to strong quantum fluctuation (Qi < 1), we can 
roughly estimate that the criterion of superfluid to insu- 
lator transition occurs when Mj unc P(Q < 1) is of order 
of one, where Afj unc is number of total junctions. From 
the data shown in Fig. H5{a), we can see that this value 
is about 9 x 0.05 = 0.45, for iV avc - log (2 x 10 5 ) = 5.3. 
Therefore we find that the quantum fluctuations of a ten- 
fragment condensate with iV avc ~ 2 x 10 5 per well should 
become strong when the disorder strength is about (or 
slightly larger than) u®. This is consistent with the esti- 
mate (see previous subsection) via the self-trapping phe- 
nomena observed in a meanfield dynamics associate with 
shaking experiments shown in Fig. |8| 

Finally we note that the results shown in Fig. I^b) are 
obtained by the same effective ID interaction strength 
used in Fig. [OJa) . In the realistic experiment, however, 
<?id will be changed if one tunes the wire current and the 
bias field, B± simultaneously (in order to keep conden- 
sate at the same position). This is because the effective 
radial confinement energy will be also increased during 
such process. However, it is easy to show that one can re- 
duce the confinement energy effectively by increasing the 
off-set magnetic field, -By , simultaneously without chang- 
ing any other system parameters. Therefore the final 
system can be kept almost the same as the one before 
tuning, except the disorder field has been increased or 
decreased independently via the composed process men- 
tioned above. 



IV. FINITE TEMPERATURE EFFECTS 

To estimate the validity of our calculation in finite tem- 
perature regime, we take the following four steps: (i) first 
we have to clarify the general concept about the tem- 
perature effects in low dimensional BEC system. (We 
do not need to include any disorder potential at this 
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point.) Due to the finite size effects in the longitudinal 
direction, it is possible that a boson system becomes ID 
quantum degenerate (or say the lowest energy eigenstate 
becomes macroscopically occupied, see Ref. when 
the temperature is below the degeneracy temperature, 
Td = Nu>\\/la(2N) > where N is the total number of atoms 
and W|| is the confinement frequency in the longitudi- 
nal direction. When considering the interaction-induced 
thermal phase fluctuation effects for temperature below 
Td, Petrov et. alrk show that a true condensate, where 
both density and phase fluctuations are small within the 
finite system size, can be achieved when temperature is 
below another temperature scale, T p h — x lju /fi, where 
/i is the chemical potential. For a temperature between 
these two temperatures, i.e. T p h < T < Td, the system 
has a frozen density fluctuation but finite thermal phase 
fluctuation. The typical experimental parameters (as we 
consider in this paper in a superfluid regime), we can esti- 
mate that T d ~ 13.2 fiK, and T ph ~ 127 nK (for N = 10 6 , 
wii = 2ir x 4 Hz, and that chemical potential is about 3 
kHz for atom- wire separation being about 100 /Ltm). The 
present experiments were done at a temperature of about 
100 nK 9 . Therefore, we can safely say that the conden- 
sate we consider here (similar to the present experimental 
conditions of MIT group) is in the true condensate regime 
(or deep in the quasi-condensate regime), where thermal 
fluctuation is very small or the thermal fluctuation length 
is comparable to the whole system size. 

(ii) Secondly, we consider the presence of a strong dis- 
order potential, where we approximate the fragmented 
condensate by a single band Bose-Hubbard model and 
solve its dynamics in meanfield approximation. The first 
approximation can be easily justified in our system from 
the following two reasons: first, we can turn on the dis- 
order potential (by moving the atom cloud closer to the 
wire) adiabatically, so that the ground state remains in 
condensate without high-energy excitations. Actually, 
the effective temperature compared to the band gap be- 
comes even smaller due to the shrink of band energy, 
since no thermal reservoir is connected to the atom cloud 
trapped in the magnetic field (thermal fluctuation from 
the wire surface^ can be neglected here because of the 
larger atom- wire separation). This is very similar to the 
situation in an optical lattice. Secondly, each fragment 
itself is in the true condensate regime due to large number 
of atoms (~ 10 4 ~ 5 ) and stronger confinement provided by 
the correlated disorder potential. The critical tempera- 
ture estimated for a single fragment is about T p h ~ 1 fiK, 
well above the temperature quoted in the experiments. 
Therefore, the thermal excitation to a higher energy band 
is strongly suppressed and can be safely neglected. 

(iii) However, the second approximation, the meanfield 



approximation for the dynamical motion of condensates, 
is self-justified only when in the deep superfluid regime 
and when the temperature is lower than both the Joseph- 
son energy Kjji and the local chemical potential devia- 
tion, Sfj, = UjiNj ~ 7V°) (see Eqs. which are 
the only two energy parameters in a meanfield version 
of the Hubbard model. Using the typical parameters we 
tracked from the data of the MIT group, the above cri- 
teria can be safely fulfilled due to the large number of 
atoms per fragment. 

(iv) Finally, we apply such a single band Hubbard 
model to study the quantum phase transition by increas- 
ing the disorder potential strength (Fig. |HJ or by reducing 
the number of atoms per well (Fig. [SJ. When considering 
the quantum fluctuation effects near the transition point, 
we note that the temperature has to be below the charg- 
ing energy, Uj, instead of NjUj in the classical limit, in 
order to see the quantum effects. Such a condition, how- 
ever, may not be satisfied in our system. In other words, 
the finite temperature effects may smoothen the transi- 
tion/crossover from the superfluid phase to the insulating 
Bose glass phase, but such crossover should still be con- 
trolled by the zero temperature quantum critical point. 
Therefore, we believe our results, like Fig. [!|]and dynam- 
ical properties near the insulating phase, should still be 
qualitatively valid at finite temperature. 



V. SUMMARY 

In summary, we study in detail the disorder magnetic 
potential in an atomic waveguide (microchip) and quan- 
titatively explain the fragmentation phenomena observed 
in the experiments. Our results shows that the dis- 
order effects can be very strong even for a very small 
wire edge fluctuation, and hence not negligible in most 
of the present experimental situations. We also study 
the dynamical properties of an array of fragments (self- 
trapping, multi-pole oscillation and modulational insta- 
bility), and propose a superfluid-to- insulator crossover in 
strong disorder limit, which can be probed by a shaking 
experiment within the experimentally accessible regime. 



VI. ACKNOWLEDGEMENT 

The authors thank useful discussion with E. Demler, 
A.E. Lcanhardt, S. Kraft, L.D. Lukin, D.E. Pritchard, 
G. Rafel, E. Altman, and W. Hofstetter. We especially 
acknowledge the critical discussion with A.E. Leanhardt 
for the experimental details. 



1 W. Ketterle and N.J. van Druten, Phys. Rev. A 54, 656 
(1996). 

2 D.S. Petrov, et. at, Phys. Rev. Lett. 85, 3745 (2000); ibid. 



87, 050404 (2001); U. Al Khawaja, J.O. Andersen, N.P. 
Proukakis and H.T.C. Stoof, Phys. Rev. Lett. 88, 070407 
(2002). 



15 



3 J.O. Andersen, U. Al Khawaja, and H.T.C. Stoof, Phys. 
Rev. Lett. 88, 070407 (2002). 

4 W. Hansel, et. al, Nature 413 498 (2001); J.H. Thywissen, 
et. al., Eur. Phys. J. D. 7, 361 (1999); R. Folman, et. al., 
Adv. At. Mol. Opt. Phys. 48, 263 (2002). 

5 A.E. Leanhardt, et. al, Phys. Rev. Lett. 89, 040401 

(2002) . 

6 D. Cassettari, et al., quant-ph/0003135 T. Schumm, et al., 
Nature Phys. 1, 57 (2005). 

7 K.K. Das, M.D. Girardeau, and E.M. Wright, Phys. Rev. 
Lett. 89, 170404 (2002); K.K. Das, G.J. Lapeyre, and E.M. 
Wright, Phys. Rev. A 65, 63603 (2002). 

8 L. Tonks, Phys. Rev. 50, 955 (1936); V. Dunjko, V. Lorent, 
and M. Olshanii, Phys. Rev. Lett. 86, 5413 (2001); M.D. 
Girardeau and E.M. Wright, Phys. Rev. Lett. 87, 210401 
(2001). 

9 A.E. Leanhardt, private communication. 

10 A.E. Leanhardt, et. al., Phys. Rev. Lett. 90, 100404 

(2003) . 

11 J. Fortagh, et. al, Phys. Rev. A 66 041604 (2002). 

12 S. Kraft, et. al, J. Phys. B, 35, L469 (2002). 

13 M.P.A. Jones, C.J. Vale, D. Sahagun, B.V. Hall, C.C. 
Eberlein, B.E. Sauer, K. Furusawa, D. Richardson, E. A. 
Hinds, cond-mat/0308434 (unpublished). 

14 T. Paul, P. Leboeuf, N. Pavloff, K. Richter, and P. 
Schlagheck, Phys. Rev. A 72, 063621 (2005). 

15 K. Huang and H. F. Meng, Phys. Rev. Lett. 69, 644 (1992). 

16 S. Giorgini, L. Pitaevskii, and S. Stringari, Phys. Rev. B 
49, 12938 (1994). 

17 A.V. Lopatin and V.M. Vinokur, PRL 88, 235503 (2002). 

18 U. C. Tauber and D. R. Nelson, Phys. Rep. 289, 157 
(1997). 

19 T. Giamarchi and H. Schulz, Phys. Rev. B 37, 325 (1988). 

20 M.P.A. Fisher, et. al. Phys. Rev. B, 40, 546 (1989). 

21 J.D. Reppy, J. Low. Temp. Phys. 87, 205 (1992); P.A. 
Crowell, et. al, Phys. Rev. B 51, 12721 (1995). 

31 B. Damski, J. Zakrzewski, L. Santos, P. Zoller, and M. 
Lewenstein, Phys. Rev. Lett. 91, 080403 (2003). 

23 J. E. Lye, L. Fallani, M. Modugno, D. S. Wiersma, C. Fort, 
and M. Inguscio, Phys. Rev. Lett. 95, 070401 (2005); T. 
Schulte, S. Drenkelforth, J. Kruse, W. Ertmer, J. Arlt, K. 
Sacha, J. Zakrzewski, and M. Lewenstein, Phys. Rev. Lett. 
95, 170411 (2005). 

24 L. Fallani, J. E. Lye, V. Guarrera, C. Fort, and M. Ingus- 
cio, cond-mat/0603655 

25 D.-W. Wang, M.D. Lukin, and E. Demler, Phys. Rev. Lett. 
92, 076802 (2004). 

26 C. Hankel and S. Gardiner in |cond-mat/0212415| C. Han- 
kel, S. Potting, and M. Wil kens, Appl. Phys. B 69, 379 
(1999). M.P.A. Jones, et. o;., |cond-m at/03010T8| 

27 J. Esteve, C. Aussibal, T. Schumm, C Figl, D. Mailly, I. 
Bouchoule, C. Westbrook, and A. Aspect, Phys. Rev. A, 
70, 043629 (2004). 

28 T. Schumm, J. Esteve, C. Aussibal, C. Figl, J.-B. Trebbia, 



H. Nguyen, D. Mailly, I. Bouchoule, C. Westbrook, and A. 
Aspect, Euro. Phys. J D 32, 171 (2005). 

29 C.J. Penthick, and H. Smith, Bose-Einstein Condensation 
in Dilute Gases (Cambridge, New York, 2002). 

30 We note that in principle, one should solve the 3D Gross- 
Pitaevskii equation (GPE) for the enlongated trap poten- 
tial, Vo(r) = ^m(bJ\r\ + lo^z 2 ), in the presence of dis- 
order potential, 5U(z) to obtain the 3D condensate den- 
sity profile, and then integrate out the transverse dimen- 
sions to get the one- dimensional density profile to compare 
with the experimental data. However, in above calculation 
we approximate such 3D calculation by the ID meanfield 
Gross-Pitaevskii equation for simplicity, where the trans- 
verse mode is assumed to be in the gaussian lowest ground 
state and the whole transverse confinement effects are ab- 
sorbed into the effective scattering amplitude, gm- Such 
approximation is valid if anio <C 1, where nm is average 
ID density. In the soldium system of MIT group we con- 
sider above, anm ~ 1.5 (for d = 100 (jm) is just above 
order of one. Therefore only few (not many) transverse 
modes are occupied, which should not be well-described 
either by 3D meanfield (good for aniu 3> 1) or by ID 
meanfield (good for anm <C 1). But since we are interested 
in strong confinement limit in all of our systems (but not 
so strong to Tonks gas limit), we think the ID meanfield 
approximation is still an appropraite approximation. 

31 D. Jaksch, et. al, Phys. Rev. Lett. 81, 3108 (1998). 

32 E. Altman, Y. Kafri, A. Polkovnikov, G. Refael, Phys. Rev. 
Lett. 93, 150402 (2004). 

33 See for example, M. Tinkham, Introduction to Supercon- 
ductivity (McGraw Hill, New York, 1996). 

34 A. Polkovnikov, S. Sachdev, and S.M. Girvin, Phys. Rev. 
A 66, 53607 (2002). 

35 A. Smerzi, et. al. Phys. Rev. Lett. 89, 170402 (2002); F.S. 
Cataliotti, et. al, New J. Phy. 5, 71.1 (2003). 

36 F.S. Cataliotti, Science 293, 843 (2001). 

37 papers for superfluid-to-Mott insulator transition. 

38 A. Polkovnikov and D.-W. Wang, Phys. Rev. Lett. 93, 
070401 (2004). 

39 T. Giamarchi, P. Le Doussal and E. Orignac, Phys. Rev. 
B 64, 245119 (2001). 

40 Y. Castin and R. Dum, Phys. Rev. Lett. 79, 3553 (1997). 

41 A. Sinatra, C. Lobo, and Y. Castin, Phys. Rev. Lett. 87, 
210404 (1997); A. Polkovnikov, Phys. Rev. A, 68, 033609 
(2003); A. Polkovnikov, Phys. Rev. A, 68, 053604 (2003). 

42 E. Altman, and A. Auerbach, Phys. Rev. Lett. 89, 250404 
(2002). 

43 S. Sachdev. Quantum Phase Transitions (Cambridge, New 
Tork, 1999). 

44 Note that in the thermodynamic limit, various regimes 
may be distinguished only for a bounded disorder, while in 
a finite system there is no true phase transition — however, 
one may still observe the sharp crossovers. 



